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Abstract 

The structure of the renormalization-group flows in a model with three quartic 
coupling constants is studied within the e-expansion method up to three-loop or- 
der. Twofold degeneracy of the eigenvalue exponents for the three-dimensionally 
stable fixed point is observed and the possibility for powers in -^/e to appear in 
the series is investigated. Reliability and effectiveness of the e-expansion method 
for the given model is discussed. 



PACS numbers: 64.60.Ak, 64.60.Fr, 75.40. Cx, 75.50.Ee 



1 Introduction 



In the this paper we consider the model given by the effective Landau- Ginzburg- Wilson 
Hamiltonian with three quartic coupling constants: 



Here ip a is a real vector order parameter field in D = 4 — 2e dimensions, a = 1,2, 
. . . , 2N, (3 = 1, 2, . . . , N. The squared "bare mass" wig is a linear measure of the 
temperature, and uo, v , and zo are the bare coupling constants. Field <p a can be 
thought of as possessing two sets of components, even and odd, each may be considered 
as an iV-component real vector. 

The Hamiltonian (P governs critical thermodynamics in a number of interesting 
physical systems. For example, when N — 2 it describes the structural phase transition 
in Nb0 2 crystal and antiferromagnetic phase transitions in TbAu2 and DyC 2 for v — z. 
Another physically interesting case iV = 3 is relevant to the antiferromagnetic phase 
transition in K 2 lrCl§ crystal and, for v — z, to those in TbD 2 and Nd |2|. The 
detailed analysis of these systems along with the Landau phenomenological theory can 
be found in Refs. |], |3|, || with references to experimental works therein. 

For the first time the renormalization-group (RG) analysis of the model (|l[) was 
carried out to the second order of e expansion by Mukamel and Krinsky in Refs. 0, 

On this ground, it was shown that the 2iV-component real anisotropic model (]l|) 
possesses a unique (three-dimensionally) stable fixed point for each N > 2. On the 
other hand, the critical behavior of this model was studied within the field-theoretical 
RG approach in three dimensions on the base of two- and three-loop approximations 
[§, §j. There were obtained expansions for (3 functions and critical exponents for 
arbitrary N. Using the generalized Pade-Borel transformation, the coordinates of the 
fixed points were found. It was shown that a stable fixed point did exist in the three- 
dimensional RG flow diagram for iV > 2. 

Assuming v = z, model (|l|) formally turns into that with generalized cubic anisotropy 
and a complex order parameter field. The latter is a specific case (m = 2) of the 
well-known mn-component model. The critical thermodynamics of this model was in- 
vestigated in detail in Refs. 0, || ||. Direct three-loop calculations for the case m = 2, 
n > 2 predict stability of the mixed fixed point, the analog of the stable tetragonal 
fixed point of model ([!]). 

In the meantime, there are general arguments, not relying upon perturbation theory, 
in favor of that a unique stable fixed point should not be in the physical space although 
its existence is not forbidden at D > 3 fllPR . The same considerations lead to the 
conclusion that the only three-dimensionally stable fixed point may be the Bose one 
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and it is that point which governs the critical thermodynamics in the phase transitions 
mentioned. The point is that when v = z, model ([!]) describes N interacting Bose 



systems. As was shown by Sak the interaction term can be represented as the 
product of energy operators of various two-component subsystems. It was also found 
that one of the eigenvalue exponents characterizing the evolution of this term under the 
renormalization group in the neighborhood of the Bose fixed point is proportional to the 
specific-heat exponent a. Since a is believed to be negative at this point, as confirmed 
by highly precise up-to-date experiments with liquid helium [12] including those in 
outer space [jl3| and the high- loop RG computations carried out for the simple 0{n)- 



symmetric model in three dimensions [14, 15], the interaction is irrelevant. Therefore, 



the Bose fixed point should be stable in three dimensions. 

However, the RG approach, when directly applied to model (^]) and to the relative 
mn-component model, has not yet confirmed this conclusion. On the contrary, all cal- 
culations performed up to now indicate the existence of a unique stable fixed point in 
the physical space, while the Bose point appears to be three-dimensionally unstable 
[1-3,5-9]. This may be a consequence of the rather crude approximations used, and 
the higher order being taken into account the closer the perturbative results could be 
to the precise results. So, the aim of the paper is to investigate the critical behavior 
of the three coupling constants model (|l|) in the next, three-loop, order in e and ver- 
ify compatibility of the predictions given by the e-expansion method with the other 
techniques. 

The main result of our study is that the unique fixed point, rather than the Bose one, 
turns out to be three-dimensionally stable within the given approximation. Calculation 
of the eigenvalue exponents of this point is a nontrivial task due to their degeneracy in 
the one-loop approximation. The analysis of the problem fulfilled in this paper shows 
that such a degeneracy results in substantial reduction of the information obtained 
from high-loop approximations. 



2 Three-loop (3 functions, fixed points, and stability 

The character of the critical asymptotics and the flow diagram structure is known to 
be determined by the RG equations for quartic coupling constants. We calculate the 
perturbative expansions for the /3-functions for arbitrary iV within massless theory 
using dimensional regularization and the minimal subtraction scheme. The three-loop 
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results follows: 
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(2) 



16(JV+4) 

12mW(864C(3) + AN + 1129) + 4wV(1440C(3) + 47iV + 1796) + 

18m; 2 z(192C(3) + 391) + 72m;2 2 (96C(3) + 103) + 2u2 3 (960C(3) + 1517) + 

1512v 3 z + 36i; 2 z 2 (48C(3) + 35) + 72u2 3 (16C(3) + 25) + 4z 4 (48C(3) + 91) 

where ( is the Riemann ( function: £(3) = 1.20206. The model under consideration 
is known to have eight fixed points @, |5|. Below we write out the coordinates of the 
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most interesting II-tetragonal fixed point only. 



u, 



N+A e + 7CTw( 70iV2 - 205iV + 139 ) r2 



(5AT-4) ' (4-5V) 3 

(iSSfC(3)(64AT 3 - 188N 2 + IhlN - 23) + 

4(4-5^)5 (6370iV 4 + 24149iV 3 - 144719iV 2 + 197208iV - 83256))e 3 , 

jg^(N - 2)e + ^±^(30^ + 25N 2 - 217iV + 166)e 2 - 

(g^k(3)(8iV 4 + 16iV 3 - 88iV 2 + 75N - 9) - i( ^_(1030iV 5 + (3) 

2751iV 4 + 46033iV 3 - 207590iV 2 + 267336iV - 109808))e 3 , 

$±(N - 2)e + ^±^(30N 3 + 257V 2 - 217iV + 166)e 2 - 

(§Sk(3)(8iV 4 + 16iV 3 - 88iV 2 + 75N - 9) - i( ^ F (1030iV 5 + 



2751iV 4 + 46033iV 3 - 207590iV 2 + 267336iV - 109808 



6 3 . 



In order to determine the character of stability of this point we should calculate the 
stability matrix eigenvalues A's. It is convenient, rather, to deal with the quantity 
y = - being a root of the reduced characteristic polynomial hereafter denoted P(y, e): 

- y 3 ( e ) + a(e)y 2 (e)-6(eMe)+c(e) = 0. (4) 

The coefficients a(e), 6(e), and c(e) are the formal series 

a(e) = a + a x e + a 2 e 2 + . . . , 

6(e) = 6 + 6ie + 6 2 e 2 + ..., (5) 
c(e) = c + cie + c 2 e 2 + 

A solution y(e) to Eq. (^) is obtained by consecutive calculating coefficients of the 
series y(e) = J2Vk£ k in corresponding orders in e. Normally, the polynomial P(y, 0) has 
three different roots y , and the derivative d y P(y , 0) does not vanish. It implies that 
the coefficients y^ are determined in the k-th order in e. This customary scheme does 
not work when the polynomial P(y, 0) has multiple roots. It takes place for the II- 
tetragonal fixed point of model (|l|) since the one-loop approximation, as will be shown 
below (see formula (||)), yields two equal values of y$. To treat this situation properly we 
have to thoroughly analyze the problem of expanding such twofold degenerate solutions. 

If coefficients a(e), 6(e), and c(e) themselves had been polynomials, the solution y(e) 
would have belonged to the class of so-called algebraic functions. Such a function is 
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analytical on the complex plane, except for a finite set of isolated points, where it has 
either poles or branchings of a finite order. The poles are associated with the zeros of 
the highest coefficient in P(y,e), while the branching points are associated with those 
values of e where P(y, e) has multiple roots. With twofold degeneracy of yo, branching 
of order 2 is possible at the point e = 0. But a two-valued function cannot be expanded 
in integer powers of e. Instead, it should be represented by a Puiseux series in powers 
of y/e |16| . Similar speculations are valid when the dependence of P(y, e) on e is of the 



formal series type since every coefficient of the expansion y(e) is determined by a finite 
number of terms in a(e), 6(e), and c(e). Therefore, whatever the perturbative order is 
chosen, there is an algebraic function coinciding with y(e) modulo higher terms. So, 
let us formulate the anzatz for y(e) as 
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= Vo + y^ 2 + yie + y^ +.... (6) 

We shall show that a solution to Eq. @ in the case of twofold degenerate roots at e = 
does exist in form (^]). We are also interested when noninteger powers in e do appear in 
and when they drop from y(e). The answer will be given by a theorem before which we 
introduce some notations. Let Z + /2 be the set of non-negative half-integer numbers. 
We assume that the infinite point oo also belongs to Z + /2. It is convenient to think of 
it as of an integer number. Denote [a, b] the interval in Z + /2 with boundaries a and b, 
i.e., the set of points / G Z + /2 satisfying inequality a < I < b. To distinguish intervals 
without one or two boundary points we use parentheses instead of square brackets. 

Theorem. In the case of twofold degenerate roots in the one-loop approximation, 
the solutions y(e) to Eq. are represented by series in powers of \fe. There is an 
alternative: either P(y, e) has two equal roots in every order of the perturbation theory 
or the solution y(e) splits at a finite step l s . Noninteger powers of e contribute to the 
expansion y(e) if and only if l s is a noninteger number. 

Let us define polynomials Ai(y ,yi, ...,yi), I G Z + /2, from the expansion 

u y iez+/2 

Another way to introduce them is as follows. Consider the e expansion of the reduced 
characteristic equation (ERCE) obtained by substitution of Eqs. (^) and (||) into 
Eq. (f|). At a sufficiently high order the coefficient before e m may be shown to be 
represented as a sum 

A (yo)y m + Ai(y ,yi)y m i + ... + Ai(y ,yi, ...,yi)y m -i + • • • = 0, m — l> I. 
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Further, let / = [0, l s ] be the maximum interval containing such that for I G / the 
coefficient yi is found from the order 21 of ERCE. It implies, in particular, that the 
partial solution [y]i = (y ,yi, ...y{) to Eq. (^) exists for all I G [0,/ s ). The upper 
boundary l s may be either finite or infinite. First let us prove that the polynomials 
Ai(yo, yi , yi) for / G [0, l s ) turn zero upon substitution of the partial solution [y]i into 
them. This is the case, at least, for I = 0. Supposing it for all I such that < / < m, 
where m is some half- integer number strictly below l s , we have 

[A (yo)y 2m +i + Ai(y ,yi)y 2m + ...] + A m (y ,yi, ...,y m )y m+ i + . . . = 

in the order 2m + \ of ERCE. The expression within the square brackets vanishes 
due to the assumption made, while the rightmost dots stand for the terms depending 
on yi with I less than m + |. Thus, if A m (y ,yi,...,y m ) ^ 0, the coefficient y m+ i is 
determined from the order 2m + \ ^ 2m + 1, in contradiction to m + \ < l s and 
m + \e I. The consequence of the fact just stated is that all yi with I G [0, l s ) are 
found from a quadratic equation 

QM = ^P(y ,0)yf + Li(y ,yi,... ,y t _i)yi + Ri{y ,yi,.. .,3/ { _i) = 0. 

The highest coefficient ^dyP(yo, 0) is nonzero as only two of the three roots yo coincide. 
Now we can give one more characteristic to the polynomial Ai(y ,yi, ...,yi) for / G /. 
Namely, it is the derivative of Qi(yi) with respect to yi. It follows directly from here 
that the solution y(e) does not split at the orders / G [0,/ s ), that is every polynomial 
Qi(yi) has equal roots. Another implication is that all noninteger / G [0, l s ) give y\ = 0. 
Indeed, every term of the linear coefficient Li(y , yi, ...,yi_i) of Qiiyi) depends on some 
variable yt with noninteger numbers k < I (because none of the series a(e), 6(e), and 
c(e) involves noninteger powers in e) and vanishes once they turn zero. Assuming 
recurrently y^ = for noninteger k < I, we find yi obeying the equation Qi(yi) = 
^dyP(yo,0)yf+Ri(yi,...,y l _i) = 0. Since Qi(yi) has two equal roots, both are zero. We 
shall show now that with finite / s , the solution y(e) splits at the order l s or, equivalently, 
that Ai a (y ,yi, ...,yi s ) ^ 0. Supposing the opposite, consider the coefficient of ERCE 
in the noninteger order 21 s + | in e: 

[Ao(yo)y 2 i a+ i + Ai{y , yi)y 2 i a + ■■■ + A la (y , y±, yi a )y ls+ i] + . . . = 0. 

The expression within the square brackets vanishes as well as the terms depicted by 
the dots on the right. Those terms contain y\ with noninteger numbers I from the 
interval /. As was stated above, such yi are equal to zero, hence in the order 21 s + | 
Eq. (|) holds identically. The next order 2l s + 1 gives a quadratic equation for y is+ i 
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with the highest coefficient |<9 2 P(j/o, 0) 7^ 0. It means that l s + \ belongs to the 
interval / that contradicts the assumption about l s being its upper boundary. So, 
Ai s (yo, yi, yi s ) 7^ and for arbitrary m 6 Z + /2 the quantity yi s + m is determined 
in the order 2l s + m in e from a linear equation in which Ai s (y , yi , yi s ) is the 
highest coefficient. Thus we have proved the existence of solutions to the reduced 
characteristic equation in the form (^). We have yet to verify disappearance of 
noninteger powers of e from expansion y(e) provided l s is an integer number. As was 
shown, the coefficient yi s + m is determined in the order 2l s + m from a linear equation 
of the form Ai s (y , yi, yi s )yi s+m + B ls+m (y , yi, % s+m _|) = 0. Supposing l s + m 
and therefore m to be noninteger, we see that each term in Bi g+m depends on some 
yi with noninteger number I < l s + m (because none of a(e), 6(e), and c(e) contains 
powers of i/e). It has already been shown that yi = for noninteger I from the interval 
/. Assuming recursively yi — for I < l s + m, we have Bi 3+m (y , yi, y m _±) = and 
therefore yi = for all noninteger Vs. 

Let us apply the theorem proved to the problem of calculating eigenvalue exponents 
for the II-tetragonal fixed point. The explicit form of the coefficients in Eq. ([5]) reads 

Ql = |_ (27(W 3 - 1129N 2 + 1591 N - 736), 

fl 2 = 2 (5iv-4)5 (48C(3)(5A^ - 4)(144A^ 4 - 72(W 3 + 1289N 2 - 9A7N + 230) 
+ 10030iV 5 - 104229N 4 + 429747N 3 - 804632N 2 + 69162CW - 222720), 

b o = (^(iV-2)(lliV-10), ( ? ) 
b x = (5j f_ 2 4 )4 (510N 4 - 3157N 3 + QQ15N 2 - 5832N + 1868), 

b 2 = 1 (48C(3)(5A^ - 4)(272A^ 5 - 1824N 4 + 4455A^ 3 - 5095A^ 2 
+ 2754N - 558) + 25890N 6 - 338437N 5 + 1547050N 4 - 3437182N 3 
+ 4044203N 2 - 2430752N + 589412) 

c o = (^(iV-2) 2 , 

Cl = ^_ 4 (N - 2)(150A^ 3 - 809N 2 + 1229N - 566), 

c 2 = 2 (5JV-4)« ( 48 C(3)(A r - 2)(5iV - 4)(80A^ 4 - 464A^ 3 + 865AT 2 

- 641N + 164) + 13950N 6 - 184745N 5 + 887705N 4 - 2072060N 3 

+ 2541094N 2 - 1575640N + 389512). 

Substituting this into Eq. (||) and setting e = (one-loop approximation) we find yo 
to be twofold degenerate: 

2-N 



(1) 1 (2) (3) 

Uo = -1, 2/0=2/0 



5JV-4' 



The solution y(e) expanding simple root is calculated in the conventional way. 
Our further consideration concerns the multiple root only. The first appearance of the 
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coefficient yi occurs in the order | in e of ERCE, with multiplier A (y ) = d y P(yo, 0): 

yiA (y ) = 0. 

Due to the degeneracy of yo, A (y ) = and yi cannot be actually determined from 
the order \. To find yi we must solve the quadratic equation in the order 1 of ERCE: 

* 2 

Qi(y±) = ld 2 P(yo,0)yl + (a 1 y 2 -b 1 y + c 1 ) = 0. (9) 



The highest coefficient ^d 2 P(y ,0) = —(3y — a ) is nonzero because only two of the 
three roots yo coincide. Substitution of Eq. ([?D into Eq. (||) gives yi = for all N. 
The next order | of ERCE does not provide y\ because at this step the equation holds 
identically: 

A (yo)ya +Ai(y ,yi)y 1 + (...)yi = 0. 
Here Ai(y ,yi) = d yi Qi(yi) = —2(3y — a )yi = 0. Considering the factor before e 2 

2 2^22 2 

in ERCE, we come to the quadratic equation 

QiiVi) = ^d 2 y P(y ,0)y 2 + (2a 1 y -b 1 )y 1 + a 2 y 2 -b 2 yo + c 2 = 0, (10) 
which has the solution 

3(iV - 1)(4(W 3 - 208iV 2 + 2537V - 66) 
Vl ~ (2N- l)(5iV-4) 3 

4|(iV-l)(iV-2)(iV + 4)(5iV-4)| 

(2N- l)(5A^-4) 3 ' 1 ' 

The order | of ERCE gives rise to the relation 

[A (y )yi + Ai(y ,yi)y 2 ] + Ax(y , yi, y x )yz + (...)yi = 0, 

where the expression within the square brackets is zero and 

Ax(y , yi, yi) = -fyoVi + ^oVi + la-iVo - h- (12) 

As was shown when proving the theorem, Ai(y , y±,yi) is the derivative of the quadratic 
polynomial Qi(yi) with respect to y±, formula (0). Hence, Ai(y , yi, yi) vanishes if 
and only if Qi(yi) has two equal roots. It occurs when N — 1 and iV = 2 in formula 
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Thus, for N ^ 1, 2 the solution y(e) splits at the integer order in e and, according 
to the theorem, would not contain noninteger powers of e. We cannot make such 
a statement for the physically significant case N = 2 because the coefficient yz is 
determined in the third order in e of ERCE involving four-loop contributions, li y3 
comes out to be nonzero, the expansion y(e) will contain noninteger powers. Otherwise, 
yz = and only the five-loop approximation, as follows from the theorem, will allow 
us to calculate the next coefficient y 2 in Eq. (^). 

To sum up, the three-loop eigenvalue exponents for the Il-tetragonal fixed point 
read 



A l = ~ e + (5iV-4)W-l) (60iV3 ~ 160N2 + 181N ~ 85)e2 + 

2 ( 5j v-4)4(i- 2 v)3 ( 4 8C(3)(2iV - l) 2 (5iV - 4)(32iV 4 - 128iV 3 + 212iV 2 - 
153iV + 33) + 20560iV 7 - 165328iV 6 + 644392iV 5 - 1406864iV 4 + 
1756745iV 3 - 1224341iV 2 + 433704iV - 59052)e 3 , 

a 2 = M e + (sjv-ir^-D ( 4 s 9< N - ^i 5 ^ 3 + 6iy2 ~ A8N + 32 i- 

3(40A^ 3 - 208A^ 2 + 253iV - 66)) e 2 , 
As = M~ 4 e + (^jgr^.!) (4 sgn(N - l)|5iV 3 + 6iV 2 - 48iV + 32| + 
3(40iV 3 - 208A^ 2 + 253iV - 66)" 



(13) 



For the physically interesting cases N = 2 and N = 3 they are presented in Table I. The 
eigenvalue exponents for the Bose fixed point, the calculation of which is an easy task, 
are also written out for comparison. It follows from the table that the Il-tetragonal 
fixed point is absolutely stable in three dimensions (3D), in contrast to the Bose one. 
Obviously, simple resummation procedures, such as the Pade and Pade-Borel methods, 
being applied to A's, do not change the picture. 

Let us now calculate the critical dimensionality iV c of the order parameter. Its e 
expansion is found from condition v c = z c = imposed on the right-hand side of Eq. 

©: 

N c = 2-26+^(6C(3)-l) e 2 + 0(6 3 ). (14) 
6 

The critical dimensionality separates two different regimes of critical behavior of the 
model. For N > N c the Il-tetragonal rather than the Bose fixed point is three- 
dimensionally stable in 3D. At N = N c they interchange their stability so that for 
N < N c the stable fixed point is the Bose one. 
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Table 1: Three-loop eigenvalue exponents for the Bose and the II-tetragonal fixed 
points (e = ~ corresponds to the physical space) 



Type 
of fixed 
point 


N =2 


AT = 3 


Bose 


\ _ l r 14 r 2 | 768C(3)-311 3 
Au 5 fc 25 fc " r 625 

\ _ 1, | 2,2 768<(3)+29 3 
Al - 5 e +5 C 625 6 

x r | 6 r 2 384C(3)+257 3 
A 2 - e +5 C 125 6 


II-tet- 
ragonal 


\ _ , | 13 r 2 84<(3)+65 3 

A 2 - -36 

A3 — —36 


\ _ , I 58 r 2 3(123600C(3)+71621) 3 
A l " fc ~ r 55 C 166375 

\ — 1 c 2 2 
A 3 - — TT e + 605 6 



Because the series (|14D is alternating, it can be resummed by means of the Pade- 
Borel method, the result being 

26 2 46 3 / 2b\^(2b\ 

N c = a + ^exp )Ei — . 15 

c c z e \ ce J \ce J 

Here a, b, and c are the coefficients before e°, e 1 and e 2 in Eq. (0), respectively, and 
Ei(x) is the exponential integral. Setting e = |, from Eq. (|15|) we obtain the value of 
critical dimensionality 

A/" c = 1.50. (16) 

Since N c lies below two, within the given approximation the critical behavior of model 
([l]) for N = 2 and A^ = 3 must be governed by the II-tetragonal fixed point. It confirms 
the deductions given by the analysis of the eigenvalue exponents. 



3 Discussions 



In conclusion, let us briefly discuss the results of the present investigation. The struc- 
ture of the RG flows of the three quartic coupling constants model was studied within 
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the e-expansion method. The three-loop series for (3 functions of the model were 
obtained. Eigenvalue exponents for the most intriguing II-tetragonal and the Bose 
fixed points were calculated for arbitrary N. For the physically interesting values 
N = 2 and N = 3, the II-tetragonal rather than the Bose fixed point was shown to 
be three-dimensionally stable in 3D, and the critical dimensionality N c = 1.5 found 
has confirmed this conclusion. Consequently, the critical thermodynamics of the anti- 
ferromagnetic phase transitions in such substances as TbAu 2 , DyC 2 , K 2 IrCl 6 , TbD 2 , 
and Nd as well as the structural phase transition in Nb0 2 crystal should be controlled 
by the II-tetragonal fixed point. It agrees with the results by Mukamel and Krinsky 
[p], 0, |3| as well as with calculations performed within the field-theoretical RG approach 
in 3D (Refs. [||, ||) but contradicts to the non-perturbative inferences ||10|| . The dis- 



tinction of 3-D RG predictions from those of the precise theory may be regarded as 
an effect of low-order approximations. The point is that the three-loop analysis in 3D 
(Refs. [H [J) shows the two rival fixed points (Bose and II-tetragonal) to be close to 
one another, and it is natural to expect that taking into account next perturbative 
terms may change their character of stability properly. Such speculations, however, 
do not suit the e-expansion because we can judge about the closeness of points just 
indirectly, i.e. by comparing their critical exponents. The e-expansion analysis yields 
the critical exponents of the II-tetragonal fixed point to be considerably distant from 
those of the Bose fixed point [l], ||, H, |17|], and hardly one can hope that longer e series 
even resummed will bring them close to one another. All this allows us to raise the 
question whether the e-expansion method is reliable for the given model. 

The twofold degeneracy of the solutions to the characteristic equation for the II- 
tetragonal fixed point in the one-loop approximation worsens the situation. According 
to the analysis performed, the eigenvalue exponents should be represented as series in 
powers of y/e instead of e. Even if we adopt the idea that noninteger powers actually 
drop from the expansions, such a degeneracy decreases the accuracy expected within a 
given approximation. Namely, in the frame of three-loop approximation we effectively 
obtain two-loop-like pieces of the series, and to evaluate the next term (of order e 3 ) we 
must take into account the five- loop contributions [jrj|. So, computational difficulties 
grow faster than the amount of essential information one may extract from the high- 
loop approximations. This leads to the conclusion that the e-expansion method is not 
quite effective for the given model. 

Of course, the e-expansion method is perfectly reliable in the four-dimensional space 
time, i.e., in studying field systems. However, as the present investigation shows, 
one should be careful when applying it to complicated three-dimensional models of 
statistical physics, especially if insufficiently high approximations are used. 
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